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The lattice regularization of QCD provides us with the most systematic way of com- 
puting non-perturbative properties of hadrons directly from the first principles of QCD. 
The recent rapid development of parallel computers has enabled us to start realistic and 
systematic simulations with dynamical quarks. In this paper, I report on the first results 
from recent systematic studies on the lattice with dynamical quarks. 

1. INTRODUCTION 

We formulate QCD on a 4-dimensional hyper-cubic lattice with a finite lattice spacing 
a. When the lattice volume is also finite, the theory is completely finite and well-defined, 
so that supercomputers can be used to calculate non-perturbative quantities. Numer- 
ous developments and ideas in the last two decades, both in algorithms for numerical 
computations and in the computer technology, have made lattice field theory one of the 
most powerful tools in computing the non-perturbative properties of QCD. In particular, 
the development of parallel computers in the late 90's enabled us to start realistic and 
systematic simulations of QCD. 

In order to extract predictions for the real world from the numerical results obtained on 
finite lattices, the following extrapolations are required: Because the continuum physics 
is defined in the limit of large lattice volume and vanishing lattice spacing, it is necessary 
to 

1. ) extrapolate to vanishing lattice spacing, 

2. ) extrapolate to large lattice volume. 

Here, computations for the quark part contain inversions of a huge quark matrix, whose 
size is, for example, about 10 7 x 10 7 for the case of a 32 4 lattice. Because the condition 
number of the quark matrix diverges as we decrease the quark mass, a calculation with 
light quarks is quite time-consuming. Actually, even with the latest supercomputers with 
TFLOPS level computational speeds, a direct simulation with the physical u and d quarks 
is not possible. Therefore, at present, we also have to 

3. ) extrapolate to light u and d quarks, 
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Figure 1. Quenched light hadron spectrum with degenerate u and d quarks [l 



using data at around the s quark mass region. 

To obtain a reliable result for continuum physics, it is essential to have a good control 
of the systematic errors due to these extrapolations. A systematic approach is required. 
Because of the huge amount of computations required, major calculations have been 
made in the quenched approximation, in which the effects of dynamical quark loops are 
ignored. A recent extensive study of quenched QCD by the CP-PACS Collaboration [[L]], 
however, has shown a limitation of the quenched approximation. As a next logical step, 
systematic "full QCD" simulations removing the quenched approximation have started. In 
this paper, I summarize the latest status of lattice QCD simulations, choosing "dynamical 
quark effects" as the keyword. 

In Sees. |2] and 0, the latest results for the light hadron spectrum from lattice QCD are 
summarized. This provides us with the most direct test of QCD in the low-energy non- 
perturbative domain. In Sec. f|, lattice results for the mass of light quarks are discussed. 
Sec. [|is devoted to the status of the U(l) problem on the lattice. In Sec. H, lattice results 
for B meson decay constants are summarized. We find that, in all of these physics, 
dynamical quarks have significant effects. A short conclusion is given in Sec. [7]. 



2. QUENCHED LIGHT HADRON SPECTRUM 

The computation of the hadronic mass spectrum directly from the first principles of 
QCD is one of the most important tasks of lattice QCD P||. If the lattice spectrum 
reproduces the experimental data, this will provide us with the most convincing demon- 
stration of the validity of QCD as the fundamental theory for strong interactions in the 
low-energy non-perturbative domain. Furthermore, through these calculations, we can see 
how the lattice QCD describes the real physics and what is required to achieve high accu- 
racy suppressing lattice artifacts. Another motivation is a determination of fundamental 
parameters in QCD; see Sec. £|. 
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Because huge computer power is required to overcome the systematic errors discussed 
in Sec [I], major simulations have been done in the quenched approximation. With this 
approximation, the amount of computations is reduced by a factor of several hundreds. 

Recently, the CP-PACS Collaboration performed a large-scale calculation of this with 
a precision significantly better than previous studies 0. Using a dedicated massively 
parallel computer CP-PACS developed at the University of Tsukuba ||, they studied 
four lattices 32 3 x 56 to 64 3 x 112 with lattice spacing in the range a « 0.1-0.05 fm. 
The spacial lattice size was fixed to be about 3 fm with which the finite size effects are 
estimated to be maximally 0.5% in the spectrum. On each lattice, five quark masses, 
corresponding to the mass ratio of pseudo-scalar and vector mesons mps/my ~ 0.75-0.4, 
were studied. The light u and d quarks were treated as degenerate. Their final results in 
the continuum limit are summarized in Fig. 0. The light quark mass m U( i and the lattice 
scale a were fixed using the experimental values for m n and m p as inputs, while the s 
quark mass was fixed either by m K (if-input) or (0-input). See for details of the 
simulation and extrapolations. 

From Fig. [I], we see that, although the global pattern of the experimental spectrum 
is reproduced, there remain systematic discrepancies of up to about 10% (7 standard 
deviations) between quenched QCD and experiment, for both choices of the input for 
the strange quark mass. Because other systematic errors are well under control, these 
discrepancies might be due to the quenched approximation. 

3. FULL QCD SPECTRUM 

In the previous section, a limitation of the quenched approximation was made clear. 
The next logical step is to perform a "full QCD" calculation removing the quenched 
approximation, to see if these discrepancies disappear with dynamical quarks. However, 
a naive extension of quenched simulations is difficult as the computational power required 
increases several hundred fold. Recently, there was much progress in improving lattice 
actions, with which lattice artifacts due to finite lattice spacing are reduced. Therefore, 
a coarser lattice may be used for continuum extrapolations with correspondingly less 
computer time. 

The CP-PACS Collaboration adopted the combination of an RG-improved gauge action 
and a "clover" -type improved Wilson quark action, and carried out the first systematic 
investigation of full QCD. Two flavors of dynamical quarks, to be identified with degen- 
erate u and d quarks, are simulated . The heavier s quark was treated in the quenched 
approximation. 

We may ask ourselves if we can actually see an effect of dynamical quarks on the present 
lattice. The most direct signals of dynamical quarks will be decays, such as p — ► rnr, and 
the breaking of the confining string for the static quark potential. It turned out that 
these phenomena are inaccesible from the current simulation parameters (quark mass, 
lattice size, etc.) ||. We can, however, see an effect of dynamical quarks in the short 
range physics: The phenomenological string tension o ~ (440MeV) 2 from the charmonium 
spectrum, is considered to represent the physics at a scale /i ~ 400-600 MeV. Therefore, 
when we compare the quenched and full QCD at the same lattice string tension, we 
are connecting the two theories at p « 400-600 MeV. Because the running of effective 
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Figure 2. Coulomb coefficient in full (filled 
circle) and quenched QCD (open circle). 
Solid line represents an estimate expected 
for the full QCD value in the chiral limit, 
and dashed lines its error. 



Figure 3. Continuum extrapolation of vec- 
tor meson masses and rnx* in full and 
quenched QCD, using the K meson mass as 
input |J. 



coupling is different between quenched and full QCD, we should see a difference in the 
Coulomb coefficient a which represents the physics at fi ~ 1-2 GeV ||. In Fig. 0, results 
for a obtained on a 24 3 x 48 lattice are shown for full and quenched QCD at the same 
lattice string tension [||. As expected, the full QCD results are larger than the quenched 
value at mpg/my = oo, and increase as the sea quark becomes lighter. The horizontal 
line in Fig. ^| is an expected value for a at mps/my = in full QCD, estimated from the 
quenched value using the two-loop beta function with two flavors of massless quarks. The 
full QCD values for light sea quarks are consistent with this estimate. Therefore, we do 
have dynamical quark effects with the expected magnitude. 

Now, we study implications of dynamical quarks in non-perturbative physics. Figure ^| 
shows the lattice spacing dependence of hadron masses from N F = 2 full QCD, compared 
with quenched masses. Plotted together are the results of a quenched calculation using the 
same improved action as in the full QCD calculation, with which one-to-one comparisons 
between quenched and full QCD results from calculations with similar parameters are 
made possible. We see that both quenched calculations lead to universal values for hadron 
spectrum in the continuum limit, and the discrepancies from experiment in quenched 
spectrum are much reduced by two flavors of dynamical quarks. The remaining difference 
might be caused by the quenching of the s quark. 

4. QUARK MASSES 

Massses of quarks belong to the most fundamental parameters of QCD. Because quarks 
are confined in hadrons, their masses have to be indirectly inferred from hadronic ob- 
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Figure 4. Continuum extrapolation of the average u and d quark mass m U d and the s 
quark mass m s in the MS scheme at 2 GeV. m s is from the .fT-input. Filled symbols are 
for full QCD. Quenched results with the standard action and the improved action are 
shown with thin and thick open symbols, respectively. 



servables using the functional relation between these hadronic quantities and QCD pa- 
rameters. A lattice QCD determination of the hadron spectrum provides us with such 
functional relations without resorting to a phenomenological model. 

On the lattice, there exist several alternative definitions for quark mass. They should 
converge to a universal value in the continuum limit. The CP-PACS Collaboration ob- 
served that light quark masses determined from either the vector Ward identity (VWI) 
or the axial- vector Ward identity (AWI), while differing at finite lattice spacings, actu- 
ally converge to a common value in the continuum limit [p]]. See Fig. [|. However, in 
the quenched QCD, the resulting value for m s differs by about 20% between i^-input 
and 0-input, due to the discrepancy of quenched spectrum from experiment discussed in 
Sec. |. 

In Fig. |], Np = 2 full QCD values, as well as quenched values with the improved action, 
are also plotted JJJ . We again find that the two quenched simulations lead to a universal 
value in the continuum limit. Results for the light quark masses in the continuum limit 
are summarized in Table [I]. The discrepancy between the K and 0-inputs becomes much 
smaller in full QCD. The most impressive point is that the values predicted through full 
QCD are 20-30% smaller than those in the quenched QCD. In particular, the s quark 
mass in two-flavor full QCD is about 90-100 MeV, which is significantly smaller than 



Table 1 

Results for quark masses in the MS scheme at 2 GeV 

m ud (MeV) m s (MeV) (if-input) m s (MeV) (0-input) 

N f = standard 4.55(18) 115(2) 143(6) 

N f = improved 4.4(2) 110(4) 132(3) 

Nf = 2 3.5(3) 89(7) 99(5) 
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the phenomenological value ^150 MeV often used in hadron phenomenology, and almost 
saturating an estimate of the lower bound 90-100 MeV from QCD sum rules M. 

5. U(l) PROBLEM 

Another important issue in QCD is the U(l) problem, that is to clarify the mechanism 
for a large rf meson mass. Propagators of flavor non-singlet meson consist of a loop of 
valence quark propagator, while propagators of flavor singlet r( meson have an additional 
contribution with two disconnected valence quark loops. The fact that rj is much heavier 
than the corresponding non-singlet meson 7r means that the two-loop contribution should 
exactly cancel the 7r pole of the one-loop contribution, leaving the heavy rf pole. This 
phenomenon is considered to be related with the anomalous violation of the flavor singlet 
axial U(l) symmetry and topological structure of gauge field configurations in QCD. This 
is a critical test of QCD and can be answered using lattice QCD. 

Calculation of the two-loop contribution requires a large amount of computations. For 
this reason only limited results have been available. The first systematic investigation in 
full QCD including continuum extrapolations was recently done by the CP-PACS Col- 
laboration ||. In an approximation ignoring with the ss state, the flavor-singlet uu + dd 
meson mass was estimated as m uU+dd = 863 ± 86 MeV. In the real world, the uu + dd 
state mixes with the ss state to lead to r](5A7) and r/(958) mesons. Therefore, the result 
863 ± 86 MeV which is slightly smaller than the experimental value 958 MeV is quite 
encouraging. Further studies, including an inspection of the mixing with the ss state as 
well as the study of topological structures, are under way. 



6. B MESON DECAY CONSTANTS 

The study of the non-perturbative properties of QCD is important also to investigate 
the Weinberg-Salam theory for the electro- weak interactions [§]. In order to determine 
CKM parameters for quarks, we need to know QCD corrections (hadronic matrix elements, 
form factors, etc.) in weak decays and mixings. In particular, properties of heavy mesons 
consisting of a b quark are being studied intensively by the B factories. The precise 
determination of the decay constants and B-parameters for B(= Bj) and B s mesons will 
lead to a strong constraint on the value of CKM parameters. The lattice calculation of 
these quantities are beginning to contribute in reducing the ambiguities in these studies 

0- 

Major calculations have been done in the quenched approximation. Recently, full QCD 
calculations of the decay constants fs and fs s have been made by several groups. In this 
section, I concentrate on the dynamical quark effects in these decay constants. 

On the lattice, simulation of the heavy b quark is not a trivial extension of light quark 
simulations, because ~ 4 GeV is larger than the lattice cutoff ~ 1-4 GeV to date. 
Two methods have been developed to simulate heavy quarks directly. One is based on a 
non-relativistic effective theory of QCD (NRQCD) defined through an expansion in the 
inverse heavy quark mass [|llj] . Another employs a relativistic action and reinterpret it in 



terms of a non-relativistic Hamiltonian 12 



In Fig. D, I collect recent lattice results for fs from quenched [|13|-[2(J and full QCD stud- 
ies H20H221 . In the figure, "Extrap." is a result of an extrapolation from light quark mass 
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Figure 5. Recent lattice results for the B 
meson decay constant from quenched simu- 
lations (Np = 0) and two-flavor full QCD 
simulations (Np = 2), with various methods. 



Figure 6. Constraints on the CKM param- 
eters p and fj using lattice results. Allowed 
region from all the constraints is shaded 
by thick lines. See text for details. 



region, and "NRQCD" and "relativistic" are the results of the lattice NRQCD method 
and the relativistic method. Estimated statistical and systematic errors are included. 
Although the errors are still large, we can see that both methods are approximately con- 
sistent with each other. Comparing the quenched results (Nj = 0) and the full QCD 
Nf = 2 results, the latters seem to be about 10% larger than the former values. From 
this figure, I estimate fs = 200 ± 20 MeV. Also for the case of fs s , increase of a similar 
magnitude is observed. 

In Fig. the implication of the lattice results for and fs 3 on the unitarity triangle 
of CKM parameters is summarized for p = p(l — A 2 /2) and fj = 77(1 — A 2 /2) using the 
Wolfenstein parameterization (see || for details). Regions surrounded by two arcs are the 
allowed regions from each constraints, except for the two thin dashed lines denoted as AM S 
which show the range of lower bounds for p. In this figure, the quenched lattice results 
for Bk, Bb, etc. are used in the constraints from e and |Ki&/Kf>| PB|, while two-flavor full 



QCD results for fs and fs s are used for the constraints from AMa and AM S . Allowed 
region from all the constraints is shaded by thick lines. The increase of Jb due to the 
dynamical light quarks causes a shift of the allowed region from AMd towards positive 
p, leading to a more severe constraint. The region shaded by thin lines is the allowed 
region from the conventional phenomenological estimates for hadronic matrix elements 
P,|23|. Hence, lattice results are contributing to reduce errors in the determination of 
CKM parameters. 



7. CONCLUSIONS 



Due to the large computer power required, major lattice simulations have been made 
in the quenched approximation, in which the dynamical sea quark effects are ignored. 
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Recently, it became possible to start realistic and systematic simulations of QCD with 
dynamical quarks, through the development of parallel computers. From the first simu- 
lations, it turned out that dynamical quark effects are quite important. The effect of two 
flavors of dynamical light u and d quarks is as large as about 20% in the values of light 
quark mass and about 10% in B meson decay constants. Both of the shifts has significant 
implications to phenomenological studies of the standard model. 

I thank the members of the CP-PACS Collaboration ]24|] for discussions. In particular, 
I am grateful to Hugh Shanahan for useful comments and preparation of several figures. 
I also thank Tony Thomas and other participants of the Conference for stimulative dis- 
cussions at Adelaide. This paper is in part supported by the Grants-in-Aid of Ministry of 
Education, Science and Culture (No. 10640248) and JSPS Research for Future Program. 
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